Spin textures, screening and excitations in dirty quantum Hall ferromagnets 
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We study quantum Hall ferromagnets in the presence of a random electrostatic impurity potential. 
Describing these systems with a classical non-linear sigma model and using analytical estimates 
supported by results from numerical simulations, we examine the nature of the ground state as a 
function of disorder strength, A, and deviation. Si/, of the average Landau level filling factor from 
unity. Screening of an impurity potential requires distortions of the spin configuration, and in the 
absence of Zeeman coupling there is a disorder-driven, zero-temperature phase transition from a 
ferromagnet at small A and {Svl to a spin glass at larger A or ISv]. We examine ground-state 
response functions and excitations. 
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Quantum Hall ferromagnets (QHFMs) are interest- 
ing especially as systems in which spin configurations 
and charge density are closely linked At small Zee- 
man energy and for Landau level filling factor ly close 
to unity, this link has the celebrated consequence that 
the charged quasiparticles with the lowest energy are not 
single electrons but skyrmions These bound states 
of a minority-spin electron with one or many spin waves 
may be viewed classically as topological excitations of 
an ordered ferromagnet: in this description, the devi- 
ation of local charge density from that of a filled and 
ferromagnctically polarised Landau level is proportional 
to the topological density Q of the spin configuration. 
For a clean QHFM with sufhciently small Zeeman en- 
ergy, skyrmions or antiskyrmions are introduced at zero 
temperature on varying the average filling factor from 
J/ = 1 to larger or smaller values and, at non-zero tem- 
perature, are generated thermally in pairs together with 
spin waves. For dirty QHFMs, coupling of an electro- 
static impurity potential to the charge density offers an 
additional mechanism by which spin textures may arise: 
the consequences of such a coupling are the subject of 
this paper. 

The interplay between disorder and exchange for 
QHFMs has been examined previously from several dif- 
ferent viewpoints. Fogler and Shklovskii, building on ear- 
lier discussions, have developed a mean field treatment in 
the spirit of Stoner theory, finding for odd integer v in the 
absence of Zeeman coupling a transition between ferro- 
magnetic and paramagnetic ground states with increas- 
ing disorder strength Q. They suggested that this tran- 
sition should be apparent in transport measurements, in 
which the ferromagnetic phase is characterised by spin- 
resolved Shubnikov-de Haas oscillations and the param- 
agnet by spin-unresolved oscillations. Experimentally, a 
transition of this kind is observed with decreasing mag- 
netic field strength |^ , and its sharpness suggests that its 
origin is indeed cooperative. Within such an approach, 
proposed for higher Landau levels where skyrmions are 
normally not stable, local moments are all coUinear in 



the ferromagnet and vanish in the paramagnet. By con- 
trast, near v = 1, an alternative is that a QHFM may 
respond to disorder mainly via the direction rather than 
the magnitude of its local magnetisation. Some indica- 
tions that this can happen come from calculations for the 
fully-polarised ferromagnet at weak disorder. Here, a re- 
duction in spin stiffness with increasing disorder strength 
has been interpreted by Green |^ as a precursor of a non- 
collinear phase. Moreover, even weak disorder may nucle- 
ate a dilute glass of skyrmions (over minima in the elec- 
trostatic potential) and antiskyrmions (over maxima), as 
discussed by Nederveen and Nazarov [Q. In addition, 
at intermediate disorder strength both reduced and non- 
collinear local moments emerge from a numerical solution 
of Hartrec-Fock theory for a model with Coulomb inter- 
actions and spatially uncorrelated disorder, by Sinova, 
MacDonald and Girvin More generally, the relative 
importance for dirty QHFMs of local moment reduction 
versus the formation of spin textures will depend on the 
nature of disorder. In the following we focus on textures, 
favoured by a smoothly varying impurity potential. 

To this end, consider a quantum Hall system with v 
close to unity, impurity potential V{r), electron den- 
sity p(r), and electron-electron interaction energy U(r). 
As a first step, treat screening using Thomas-Fermi the- 
ory, omitting exchange interactions and Zeeman energy. 
Within this approximation, developed by Efros for 
the comparable problem in spin-polarised Landau lev- 
els when v lies near half-integer values, the ground-state 
charge density at weak disorder is determined by the 
condition that the Hartree potential should everywhere 
match the chemical potential: /i = V{r) + J U{r — 
r')p{r')(Pr' . We are concerned with circumstances in 
which the resulting density varies smoothly on the scale 
of the magnetic length, ^b, and has fiuctuations from 
v ~ 1, Sp{r) = p{r) — (27r^|)^^, which are small: 
|(5p(r)| ^ p(r)- Restoring exchange interactions under 
these conditions results locally in maximal ferromagnetic 
polarisation of electron spins, with a direction that may 
vary in space. Denoting this direction by the three- 
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component unit vector S{r), its spatial fluctuations are 
linked to electron density via 

Spiv) = {87T)-h,jea,f3-,S"d,S^djS'' . (1) 

An exchange energy is associated with such variations: 
combining exchange (with interaction constant J), im- 
purity and Hartree contributions to the total energy, 
and choosing for simplicity a short-range interaction 
U{r — r') = Uo5{r — r'), we take as our description of a 
dirty quantum Hall ferromagnet the configurational en- 
ergy 

I (^^\VSir)f + V{r)6p{r) + ^[Sp{r)fy'r. 

(2) 

Doing so, we neglect quantum fluctuations of S{r), as 
is justified in the semiclassical limit, |VS'(r)| <^ l^^, to 
which we are already restricted. We choose for simplicity 
V{r) Gaussian distributed with amplitude A and correla- 
tion length A. (More realistic choices for U{r) and for the 
distribution of V{r) will be considered elsewhere pO[.) 

Our aim in the following is to understand the zero- 
temperature phase diagram of the model defined by 
Eq. (|^), and to characterise its ground states via their 
response functions and excitations. As an example of a 
disordered electron system, it is unusual in that there is 
an exchange gap for single-particle excitations, even if the 
ground-state spin configuration ^(r) is not ferromagneti- 
cally ordered, so that the only low-energy excitations in- 
volve collective spin modes. As an example of a ferromag- 
net with quenched disorder, the system is also unusual in 
several ways. First, the coupling to disorder leaves spin- 
rotation symmetry intact but breaks time-reversal sym- 
metry, in contrast to random magnetic fields, which break 
both symmetries, and to random exchange interactions, 
which leave both symmetries intact. Second, because of 
the same coupling, the spin system responds to applied 
electric fields: we calculate the wavevector-dependent di- 
electric susceptibility, comparing with behaviour found 
in more conventional disordered electron systems. Third, 
the link between spin and charge also endows spinwaves 
with an electric dipole moment: we calculate the spin- 
wave contribution to the optical conductivity in disor- 
dered spin states, complementing Green's results for 
the polarised ferromagnet with impurities. 

We start with a simple discussion of the phase dia- 
gram. The model is characterised by two energy scales, 
J and A, and two length scales, A and Lh = {Uq/J)^/"^. 
The last of these, which we call the Hartree length, plays 
an important role in what follows. Its significance can be 
made clear by comparing, for a skyrmion of fixed shape 
and radius i? in a clean system, the contributions to to- 
tal energy from exchange and from Hartree interactions, 
of order J and Uo/B? respectively: exchange dominates 



on lengthscales large compared with Lh, while Hartree 
interactions dominate at smaller distances. We use the 
limit Lh 3> A as a source of simplifications in analytical 
estimates, but take Lh ~ A in numerical simulations. 

Examine first the ground state for ((5/9) = as a func- 
tion of A/ J. At weak disorder (A < J) it can be pic- 
tured as a dilute glass of skyrmions and antiskyrmions, 
discussed in Ref. 0. This reflects the existence of 
a threshold [^"D, defined by |V'(r)| — An J: in most 
parts of the system (roughly, those where |V'(r)| is be- 
low threshold) ferromagnetic order is essentially unaf- 
fected by the impurity potential More precisely, 
for any S{r) one has |V5(r)p > 87r|(5p(r)| and hence 
n > /[47rJ|(5p(r)| + Vir)dp{r)]d'^r, which for V{r) be- 
low threshold everywhere implies that the ground state 
is the perfectly aligned ferromagnet with Sp{r) = 0. At 
stronger disorder (A ^ J), by contrast, screening is al- 
most perfect at short distances and 6p{r) ~ —V{r)/Uo. 
Corrections to such screening arise at and beyond the 
scale Lh, where exchange is important. We can sum- 
marise the effect of exchange by dividing the system 
into regions of area L^, finding for each such area the 
integral Q = — J V{r)/Uo (Pr, and adjusting the to- 
tal screening charge within every region to the integer 
value closest to Q. We argue that these integers are 
predominantly zero in a ferromagnetic phase, and pre- 
dominantly non-zero in a phase without ferromagnetic 
order. To see this, consider a well-ordered ferromagnetic 
phase, in which S{r) as a function of r has small fluc- 
tuations around a global direction of magnetisation. In 
this case, the net topological charge in each region has 
magnitude much less than one. Conversely, in a phase 
without such order, unit topological charge accumulates 
in a region of linear size given by the ferromagnetic cor- 
relation length. The phase boundary for the ferromagnet 
is therefore located by setting (Q^)^/^ 1. Since each 
area of size A^ contributes to Q a charge of magnitude 
X^A/Uq with random sign, 

(Q2)1/2 ^ XLhA/Uo yielding 
Ac [/o/(ALh) = J(Lh/A). For A > A^ the ground 
state has no ferromagnetic order; because within a clas- 
sical zero-temperature description spins are frozen, we 
identify this phase as a spin glass. 

Consider next the ground state at fixed A < A^, 
as a function of (Sp). Charge is introduced into the 
system for (Sp) > as skyrmions of size B. In the 
presence of Hartree energy alone, R is divergent, but 
fluctuations of an impurity potential establish an op- 
timal size, by generating random potential wells. For 
R :S> X these have an average depth A{X/R) (the case 
i? ^ A is treated in Ref. The contributions to 

the skyrmion energy from these two sources are of or- 
der Uq/R^ and — AA/i?, respectively; minimising their 
sum, we find R ^ C/o/(AA) = Lh(Ac/A). (Note that, 
since i? > Lh for A < Ac, screening as discussed in the 
previous paragraph does not alter this argument.) We ex- 
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FIG. 1. Phase diagram for Lh ~ \/2A, obtained from 
simulations. Inset: phase diagram for Lh 2> A, obtained 
from analytical arguments described in the text. 

pect ferromagnetic order to persist with increasing {Sp) 
until such skyrniions overlap, so the phase boundary lies 
at {Sp)c ~ Ljj^(A/Ac)^. In this way we arrive at the 
schematic phase diagram shown in the inset to Fig.|l|. 

In order to test these arguments, we have studied 
ground states of a lattice model numerically. In outline, 
we replace Eq. (|^) with a Heisenberg model on a square 
lattice which has nearest-neighbour ferromagnetic inter- 
actions of strength J and a topological charge defined 
on each elementary plaquette. These charges contribute 
to the total energy by a local Hartree interaction with 
coefficient Uq, and by their interaction with plaquette 
potentials, independently and uniformly distributed in 
[—A, A]. Thus A = 1 in units of the lattice spacing. We 
find low-energy states using a simulated annealing proce- 
dure similar to that of Ref . ||l^ , and can identify states 
that are equivalent under a global rotation by comparing 
total energies and also charge distributions. Using system 
sizes up to 56^ spins and annealing with 10^ Monte Carlo 
steps per spin, for a given disorder realisation in succes- 
sive runs we repeatedly reach states from a small set: we 
take lowest of these in energy to be the ground state. In 
the ground state for each disorder realisation we calcu- 
late the site-averaged magnetisation, M = [(^(r))!. We 
also determine: the ground-state susceptibility Xi from 
the response to a Zeeman field applied in the direction 
of the magnetisation; the ground-state spin stiffenss ps, 
from the energy cost of long-wavelength twists imposed 
on the spin configuration; and the wavevector-dependent 
compressibility, from the linear response of Sp(r) to a pe- 
riodic potential added to V{r). Further details of our 
methods will be described elsewhere |l0| . 

Representative results are shown in Fig. |[ Taking 
{Sp) = and Uq/J = 2, we present M, x and ps as 
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FIG. 2. Ground-state properties as a function of disorder 
strength: magnetisation M, susceptibility x, and eigenvalues 
ps of the spin stiffness tensor. 

a function of disorder strength. A/ J, using 40^ spins 
and an average over three disorder realisations. Three 
regimes are evident from the behaviour of M. For 
A/J < 1, the ground state is a fully-polarised ferromag- 
net, because we use a bounded disorder distribution. For 
1 < A/J ^ 2.5, the ground state is a partially polarised 
ferromagnet (though with an almost saturated magneti- 
sation for 1 < A/J < 1.8). And for 2.5 < A/J, the 
ground state is, we argue, a spin glass, with a small 
non-zero M arising as a finite-size effect. In support of 
this interpretation, a large peak in x indicates a phase 
transition at A/J « 2.5. Moreover, the distinction be- 
tween phases is illustrated by their spin stiffness, spec- 
ified in full by a 3 X 3 symmetric tensor. We plot the 
eigenvalues of this in Fig.||. For the fully-polarised ferro- 
magnet, rotations about the magnetisation direction do 
not alter the spin configuration, and so one eigenvalue is 
zero while the other two are degenerate For the partially- 
polarised ferromagnet, all three eigenvalues are non-zero, 
with two remaining degenerate. And at 2.5 < A/ J, mag- 
netic isotropy and finite spin stiffness in the spin glass are 
illustrated by the fact that all three eigenvalues are ap- 
proximately degenerate and non-zero. 

Repeating such calculations for {Sp) > 0, we deter- 
mine the phase diagram shown in the main panel of 
Fig. ^ which is qualitatively similar to that in the in- 
set, even though the simulations are for {Uq/JX^) = 2, 
while our theoretical arguments apply for {Uo/jX^) » 1. 
The maximum density range over which disorder may 
stabilise the ferromagnet, \{Sp\ < 10"'^ in units of charge 
per plaquette, is strikingly narrow. 

We now turn to the dielectric response of the par- 
tially polarised ferromagnet and spin glass, characterised 
at zero frequency by the wavevector-dependent dielectric 
susceptibility Xci<l) or by the compressibility ^(g), with 
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FIG. 3. Compressibility n{q) as a function of wavevector q 
for systems with J = 1 and Uq = 1 (0,+) ot Uo = 2 (0,x); 
full line; quadratic fit at small q. 

i^il) — 9^Xo(9)eo/e^. To fix definitions, take V{r) — > 
V{r) + T^i cos(q • r) in Eq. (||): the ground-state elec- 
tron density changes according to Sp{r) Sp{r) + Spi{r) 
and the disorder-averaged linear response is {6pi{r)) ~ 
—K{q)Vi cos(q • r). We can apply to K{q) the approach 
used in our discussion of the phase diagram: the Hartree 
length Ln again plays an important role, within our 
model problem in which Coulomb interactions are omit- 
ted. For q :$> ijj"^, exchange may be neglected and 
K{q) ~ Uq^. More generally, we anticipate the scaling 
form K{q) — [/(^^/(-Lh*?), with /(x) constant at large 
X. For q <C we expect that exchange dominates 
and hence that K{q) should be independent of Uq in 
this regime. In turn this implies for the scaling function 
f{x) oc at small x, so that K{q) oc (qLnyUf^^ = Jq-^ 
for q <^ ^H^- summarise, the response is that of a 
metal (^(g) constant) at large wavevectors, and that of 
an insulator (xe(9) constant) at small wavevectors p^ ]. 

We have tested these arguments using the simulation 
methods outlined above. Results for K{q) are displayed 
in Fig.^, where we compare behaviour in systems with 
J = 1 and Uq — 1 or Uq = 2, combining for each case 
data from lattices of size 40^ and 56^ in order to maximise 
wavevector resolution. The distinction is clear between 
a C7o~iiidependent K{q), quadratic in q, at small q, and a 
g-independent «((/), varying roughly as Uq^, at large q. 

Finally, we consider optical conductivity <t{lj) at fre- 
quency uj. A spin glass in the absence of Zeeman energy 
is expected to support three degenerate Goldstone 
modes. Their dispersion is linear at small frequencies 
with speed c = (ps/x)^^^- In an QHFM they have an 
electric dipole moment, which arises because they propa- 
gate in a non-collinear ground state. Their contribution 
to (t(w) is determined by a product of the mean square 
dipole moment with the density of states. We find 



where ^ is the spin correlation length. Variable-range 
hopping would mask this contribution to cr(a;) but should 
be absent if spins locally have maximal polarisation 
throughout the sample. 

In summary, we have shown how the competition be- 
tween screening and exchange determines the ground- 
state properties of QHFMs in the presence of a smoothly- 
varying impurity potential. 
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